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ANALYTIC CALCULATION OF LAUNCH VEHICLE RESPONSE TO WINDS AND 
CALCULATION OF OPTIMAL BIASED PITCH PROGRAMS 

by Janos Borsody 
Lewis Research Center 

SUMMARY 

A set of simple, closed-form equations is derived to evaluate the trajectory re- 
sponse to a wind disturbance; that is, to compute 

(1) Angle of attack caused by a wind disturbance for a trimmed vehicle (The vehicle 

is trimmed if nominal attitude is maintained in the presence of a wind disturb- 
ance. ) 

(2) Attitude bias caused by a wind disturbance if a load-relief control system is as- 

sumed (i. e. , the nominal angle-of-attack profile is maintained in the presence 
of wind) 

(3) Thrust vector deflection required to maintain nominal flight attitude in the pres- 

ence of a wind disturbance 

In deriving the closed-form solutions to the equations of motion, some simplifying 
assumptions were made. This introduces errors in the resulting angle-of-attack or 
attitude bias. However, the results obtained using the closed-form equations agree well 
with the detailed six-degree-of-freedom computer results. 

The simplified equations are also used to derive optimum biased pitch and yaw pro- 
grams based on a sample of actual wind measurements. The results presented show a 
significant improvement in launch availability (based on a constant angle- of- attack- 
times-dynamic-pressure capability of the vehicle) compared to the nominal (no bias) 
trajectory and to a biased pitch program based on a statistical (synthetic) wind profile 
considered herein. 


INTRODUCTION 

One of the many launch vehicle studies required prior to a flight is to determine the 
launch availability (the probability that the structural capability of the vehicle will not be 
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exceeded during flight) which can be expected for a vehicle with a given structural 
strength or allowable aerodynamic loading. Alternately, in preliminary design studies 
it is necessary to determine the structural strength required for a specific launch avail- 
ability. Since maximum structural loading of launch vehicles is often caused by the 
winds aloft, both of these studies require the calculation of aerodynamic loads using 
synthetic wind profiles (refs. 1 and 2) or a large sample of actual wind measurements. 

Further studies involve the possible reduction of aerodynamic loads (or increasing 
launch availability) without structural modification of the vehicle. Aerodynamic loads 
can be reduced by using a load-relief autopilot or biased pitch and yaw pregrams. The 
basic principle of a load-relief autopilot is a continuous in-flight calculation of aero- 
dynamic loads and turning the vehicle in such a way as to minimize these loads. The 
second method of reducing aerodynamic loads is by using a biased pitch program. The 
reference (nominal) pitch program is usually defined by a zero-angle- of-attack, no- 
wind trajectory simulation. The initial pitchover is selected to maximize payload or to 
satisfy an ascent aerodynamic heating constraint. To obtain a biased pitch program, 
the nominal flight attitude is adjusted to minimize the expected aerodynamic loads based 
on a sample of actual wind measurements. Biased pitch programs make use of the 
seasonal correlation of wind velocities and directions in reducing loads. This seasonal 
correlation of winds is observed in references 1 and 2. Derivation of a biased pitch 
program involves the calculation of aerodynamic loads for a large sample of actual wind 
measurements. 

Biased pitch programs were developed for the Atlas-Centaur vehicle using a detailed 
computer program and a sample of actual wind measurements. One of the procedures 
used is discussed in reference 3. This procedure gave good biased programs. However, 
it was too costly in terms of computer time to be of practical use. To reduce the com- 
puter time required, biased pitch programs for the actual flights were derived based on 
artificial wind profiles (ref. 4) which were representative of the seasonal winds. The 
procedure gave good launch availability. 

In the past, all these studies were accomplished by using detailed six-degree- of - 
freedom computer programs. This is time consuming and expensive. To reduce the 
cost and to simplify these studies, a set of simple, closed-form equations is derived in 
this report for calculating the trajectory response to a wind disturbance. Furthermore, 
these equations are used to derive biased pitch programs using a large sample of actual 
wind measurements. 

The basic equations of this report are obtained by linearizing the equations of motion 
about a reference trajectory. The simplified equations can be used to determine angle- 
of-attack change caused by a wind disturbance if the vehicle is assumed to be trimmed 
(i. e. , the nominal attitude is maintained in the presence of a wind disturbance). If a 
load-relief autopilot is assumed (i. e. , the angle- of-attack profile is minimized in the 
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presence of a wind disturbance), the equations can be used to determine the change in 
flight attitude caused by a wind* disturbance. 

Maximum loading of launch vehicles occurs in the atmosphere near the maximum 
dynamic pressure region. Loads on the vehicle are made up of axial loads and bending 
moments. For a particular vehicle, axial loads cannot be changed, since they depend 
on axial acceleration and aerodynamic drag. The bending moments, on the other hand, 
depend on the angle-of-attack history encountered during flight. Bending moments are 
proportional to the product of angle of attack and dynamic pressure. Dynamic pressure 
does not change appreciably because of the wind disturbance. Therefore, in order to 
reduce bending moments, it is necessary to reduce the angle of attack. The simplified 
equations can be used to compute the angle-of-attack change resulting from the wind 
disturbance. This angle-of-attack change, in turn, can be used to compute the attitude 
bias (using the simplified equations) necessary to maintain a nominal (no wind) angle- 
of-attack history, thus giving the best possible bias. Since the equations are simple, it 
is possible to analyze a sample of actual wind measurements in a minimum of computer 
time. 

Results obtained by using the simplified equations are presented for two vehicle 
configurations: a 260-inch solid-rocket booster with an STV-B upper stage (configura- 
tion I) and the Atlas-Centaur (AC- 15) vehicle (configuration II). The results presented 
include comparisons of angle-of-attack and attitude bias for three wind profiles using a 
detailed six-degree-of-freedom computer program and the simplified closed-form equa- 
tions. Furthermore, the simplified equations were used to derive biased pitch and yaw 
programs for the two vehicle configurations. The biased programs are based on a sam- 
ple of 100 March wind soundings. 

To evaluate the biased pitch and yaw programs, launch availability is derived by 
using the simplified equations and the sample of 100 winds (ref. 2). A comparison of 
launch availability is made for the nominal and biased pitch and yaw programs. (The 
launch availability described is based on maximum angle of attack times dynamic pres- 
sure. This assumes the vehicle has a constant angle-of-attack-times-dynamic-pressure 
capability. Of course, if this capability is given as a function of flight time, the launch 
availability can be derived by using the given vehicle capability. ) 


ASSUMPTIONS 

To obtain a closed-form solution to the equations of motion, the following assump- 
tions are made: 

(1) The vehicle is assumed to be a rigid body; no attempt is made to include or 
compensate for elastic effects. 
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(2) The inertial effect of the vehicle is assumed to be negligible; that is, the vehicle 
can be instantaneously turned to the desired attitude. 

(3) The basic equations of motion are linearized about a nominal operating point, 
which is generally chosen to be a no-wind, zer o- angle- of-attack trajectory. Since the 
analysis gives the changes from this nominal flight caused by wind disturbances, it 
requires a detailed nominal simulation to evaluate the perturbed flight. 

(4) The change in velocity from the nominal caused by the wind disturbance is as- 
sumed to be small, and its effect on the flight path angle negligible. If greater accuracy 
is required, this effect can be included. In this case, however, two first-order linear 
differential equations (rather than one as will be discussed) must be solved simultane- 
ously to obtain the change in angle of attack. 

(5) The basic vehicle parameters (thrust, weight, drag, etc. ) are assumed to re- 
main constant in some small time interval. This allows a closed-form solution (in this 
interval) of the linear first-order variable- coefficient differential equations. 

(6) The wind velocity is approximated by a straight line in some subinterval of the 
one chosen in assumption 5. Since both these intervals are arbitrary, this does not 
limit the accuracy of the results. 


ANALYSIS 

The variables and other notations used in the following discussion are defined in 
appendix A. The basic equations of motion have been derived in reference 5, and the 
equations for the flight path angle, velocity, attitude, and angle of attack are reproduced 
in appendix B for convenience and completeness. 

The basic equation of this analysis relating the attitude change A 9, change in angle 
of attack A a, and superimposed wind angle of attack « w is derived in appendix B and 
is given by 


A *0 + aj A0 = Aa + Cj A a + a» w + b^ar w (B5a) 

The superimposed wind angle of attack at any particular flight time depends on the wind 
velocity, the wind azimuth, and the vehicle's relative velocity and flight attitude. The 
change in angle of attack, on the other hand, depends on the vehicle and trajectory re- 
sponse to the wind disturbance. In the analysis, the coefficients a^, c^, and b^ are 
assumed to be constant in some time interval. This assumption is made in order that 
a closed-form solution may be obtained for this equation. 
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Attitude Bias 


In the following sections, simple closed-form solutions of equation (B5a) will be 
derived to compute the required change in flight attitude (attitude bias) that will main- 
tain a nominal angle- of- attack profile in the presence of a wind disturbance. Simplified 
equations will also be derived to compute the attitude bias for a desired change in 
nominal angle of attack. These equations are given for a particular wind or change in 
angle of attack. Their usefulness in deriving biased pitch programs is discussed in 
the section Statistical Wind Analysis and Biased Pitch Programs. 

Pitch and yaw plane attitude bias. - It is desired to maintain the nominal angle- of- 
attack profile (in general, near zero angle of attack) in the presence of winds. This 
assures that the aerodynamic loads will be minimized. Since 

Aa = a - ar n = 0 (1) 

equation (B5a) gives the desired relation of wind- induced angle of attack and the change 
in vehicle attitude 


A0 + a x A 9 = a w + b^ 
The general solution of equation (2a) is given by 


(2a) 


A9 


= + ( AB o ~ “w, o> ex P jf' a i(T)drJ - a w (|)ex P £ a^dr ] 


x [a x U) - bjd)]^ (2b) 


Assume that coefficients a^ and b^ remain constant on the interval F^; 
Fj £t.j < t < (tj + Atj)J . With this assumption equation (2b) becomes 


ie i + i - Vj+i + (A0 j ' Vi )exp( ‘ a i A V 


/•tj+Atj 

( a l - bj) jf j+1 («)exp[- a l^j+l ' ?) ] d? (2C) 
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If the superimposed wind angle of attack is assumed to be linear on the interval E-; 
E^tj < t < (t. + Atp] where Ej C F^, then 

a w, i+1 ” ®w, i^ ~ V + a v/, i 

When equation (Ble) is used, the following expressions for a? w ^ and » w i are ob- 
tained: 


O'. 


w , i 


V w sin(9 - q) n 


t=t. 


(3b) 


and 


O', 


|V W sin ( 0 " «) n 


w ’ x At. 


- cv, 


W, 1, 


-t=t i+ At. 


(3c) 


With these assumptions, the solution of equation (2c) becomes 
A0 i + 1 = «w,i+l + (A *i " a w,i )eXp( " a l At i ) " (* ‘ “) 


x * 


(®w, i " t 1 - exp( " a l At i } ] + ®w, i At ij 


(4) 


Equation (B5) can also be used to obtain the attitude bias required for a desired change 
in nominal angle -of- attack profile. To do this, the wind-induced angle of attack is as- 
sumed to be zero (o w = c* w = 0). Furthermore, if the angle- of- attack change is approx- 
imated in the same way as previously for the wind-induced angle of attack (by a straight 
line segment), the solution to equation (B5) becomes 


A0 i + 1 = Aa, i+1 + < A0 i - Aa i )ex P(-a 1 A V “ ^ “ 


A a 


Affj -) [l - exp(-a^ At.)J + Aon A^ 


(5) 
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where 


A *°i = “< Ao i + l- AOi) 

If a particular angle-of-attack history is desired in the presence of a wind disturbance, 
equations (4) and (5) can be used to compute the attitude biases required for the wind 
alone and the angle-of-attack bias alone, respectively. Since equation (B5) is a linear 
differential equation, the total attitude bias is the sum of the attitude biases due to the 
wind and to the angle-of-attack bias. The same equations can be used in computing at- 
titude bias in the yaw plane. Of course, it must be remembered that a^, bj, and Cj 
of appendix B are computed based on the yaw plane trajectory parameters. 

Final conditions . - In the previous section, equations were presented to give attitude 
bias. This bias is required to maintain a given angle-of-attack profile in the presence 
of a wind disturbance. If the nominal trajectory is assumed to be zero angle of attack, 
the bias given by these equations is the same as would be obtained by using a perfect 
load-relief autopilot. However, this bias gives large dispersions in altitude and flight 
path angle at booster separation. This can result in payload loss or excessive aerody- 
namic heating. When biased pitch programs are derived, the dispersions can be mini- 
mized by varying the initial attitude bias to maintain nominal flight path angle at booster 
separation. This is equivalent to changing the initial pitchover in the trajectory. Lin- 
earizing equation (Bid) gives 


A a - Ad •- Ay - Q! w (6) 

Since in deriving equation (4) it was assumed that the nominal angle-of-attack profile 
will be maintained (i. e. , A a = 0), equation (6) reduces to 

Ad = Ay + o- w 

To reduce the dispersions, the flight path angle is chosen to be the same as the nominal 
flight path angle at booster cutoff (i. e. , Ay^ = 0). Then 

< 7 > 


The initial attitude change is varied to satisfy the given final condition (eq. (7)). Since 
the relative velocity is large at booster engine cutoff, the superimposed wind angle of 
attack becomes small. This assures an essentially nominal burnout condition; that is, 
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the attitude, angle of attack, and flight path angle are all close to their nominal values. 
This, of course, does not preclude small deviations in altitude and downrange position. 


Simplified Angle-of-Attack Computation 

Most present-day launch vehicle autopilots are designed to maintain nominal flight 
attitude (trimmed) in the presence of aerodynamic disturbances. If the vehicle is as- 
sumed to be trimmed, and a unity autopilot is assumed, equation (B5) can be used to 
give a functional relation between the change in angle of attack and superimposed wind- 
induced angle of attack; that is, 

A'a + c 1 Aa = -(a w + b 1 a w ) (8) 


The general solution to equation (8) is given by 


l a=-a ff + (Aa 0 +a ff) Q )exp^- c 1 ( T )d T j + £ a w (Cj - b^expj^- c^rMrJdl 


( 9 ) 


If the same linear wind approximation as in the previous section on attitude bias is used 
and it is assumed that the variables Cj and bj remain constant in the interval, the so- 
lution to equation (9) is given by 



This equation can be used in both pitch and yaw planes with the appropriate coefficients. 


Thrust Vector Deflection Requirements 

Thrust vector deflection requirements caused by the wind, for a trimmed vehicle, 
may easily be computed from the linearized form of equation (Blc); that is, 
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( 11 ) 


A6 = 


*a N l 

*c T 


A a 


This equation does not involve any characteristics of a real autopilot. However, because 
most control systems are designed to maintain near-trim conditions and their response 
time is small, equation (11) gives a good approximation to the thrust vector deflection 
requirement. 


Statistical Wind Analysis and Biased Pitch Programs 

A set of biased pitch programs has been used on some of the Atlas-Centaur flights 
with significant improvement in launch availability. However, it is difficult to obtain an 
optimum biased pitch program. There have been two different procedures used to de- 
rive the necessary bias; namely, biasing the nominal (zero angle of attack) flight to a 
statistical (synthetic) wind, or running a set of real wind soundings on a detailed com- 
puter program and adjusting the attitude (more or less empirically) until the resulting 
trajectories are acceptable in terms of aerodynamic loads. The first procedure has the 
disadvantage of giving a bias which may not be the best available. The second procedure 
gives good results. However, it requires prohibitively large computer times to derive 
a biased pitch program. 

This section presents a simple, straightforward procedure for deriving a biased 
pitch program. Since structural loading is proportional to the angle of attack at a par- 
ticular flight time, these loads can be minimized by minimizing angle of attack. The 
desired angle- of-attack bias can be derived from a statistical analysis of a set of real 
wind soundings. Equation (10) gives the change in angle of attack caused by a wind along 
the trajectory. Use of this equation does not require a trajectory integration to deter- 
mine the change in angle of attack caused by the wind disturbance. It requires the inte- 
gration of a single, first-order differential equation. In fact, if the trajectory param- 
eters are assumed to be constant, a closed-form solution is available for each interval. 
By using this closed-form solution, the angle- of-attack profiles of a large sample of 
winds can be derived in a minimum of computer time. To obtain the best angle-of-attack 
bias, a statistical average of the angle-of-attack biases (for the sample of winds) is 
taken at every time point along the trajectory. Once this angle-of-attack bias is known, 
equation (5) can be used to compute the optimum attitude bias. 

Generally, there will be large angles of attack in the early part of the flight. This 
is caused by the low relative velocity of the vehicle. The large angles of attack may 
give large angle-of-attack bias. The large angle-of-attack bias, in turn, would give 
undesirably large attitude bias. To eliminate the large attitude bias, the angle-of-attack 
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bias is set equal to zero in the early portion of flight. The angle- of-attack bias can be 
set equal to zero in this flight region without adverse effects on loads since the dynamic 
pressure is very low. 


RESULTS AND DISCUSSION 

Figure 1 presents the wind profiles used in evaluating the accuracy of the simplified 
equations. The winds were chosen to give large perturbations: that is, high angles of 
attack and wind shear. The two winds given in figures 1(a) and (b) are designed to give 
high wind velocities and high wind shear near the maximum dynamic pressure region. 
They are given in terms of wind velocity and direction as a function of flight time. This 
choice of independent variable was made for convenience. A wind azimuth of 285 de- 
grees for winds 1 and 2 was selected to eliminate any wind component in the yaw plane. 
Both vehicle configurations were launched at a 105-degree azimuth. The wind in fig- 
ure 1(c) was selected from a sample of actual wind measurements (ref. 2) based on its 
high wind velocity, high wind shear, and a broad wind plateau near maximum dynamic 
pressure. The effect of wind direction is also included in the analysis of this wind. 
However, only the pitch plane results are presented. 

Figures 2 and 3 give the angle- of-attack change from the nominal obtained by using 
both the simplified procedure and the six-degree- of-freedom computer simulation for 
two vehicle configurations, respectively. The two vehicle configurations are the 260- 
inch solid-rocket booster with an SIV-B upper stage (configuration I) and the Atlas- 
Centaur (AC- 15) vehicle (configuration II). The results obtained by using the simplified 
equations for the three wind profiles show good agreement with those obtained from a 
detailed computer simulation. 

Attitude biases were also derived for the winds presented in figure 1. In deriving 
the attitude biases, the simplified closed-form equation was used, and a nominal angle- 
of-attack history is assumed in the presence of the wind. 

Figures 4 and 5 represent the attitude biases for the two vehicle configurations. 

The curves show the bias with and without iterated end conditions. Both biases were 
obtained by changing the attitude along the trajectory to cancel the wind-induced angle 
of attack. For the noniterated bias, the initial attitude bias was set equal to zero. 
However, for the iterated bias, the initial attitude bias was adjusted to achieve the nom- 
inal flight path angle at booster cutoff. These figures show a significant reduction in 
maximum bias requirement when the iterated end conditions are used. This reduction 
is caused by an increased attitude bias in the low-wind- velocity region. Because of the 
inertial and aerodynamic properties of configuration n, the large wind velocities en- 
countered in the early part of the flight with wind 3 caused the attitude to diverge in the 
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case of noniterated bias. Therefore, the data for configuration II with wind 3 are not 
presented in figure 5 and later figures. 

Since the biases are designed to maintain nominal angle of attack in the presence of 
winds, the accuracy of the simplified equations can be evaluated by simulating a detailed 
trajectory with a wind and its corresponding attitude bias. The resulting angle of attack 
is compared with the nominal; the difference in angle of attack is caused by the approxi- 
mations made in deriving the equations. The difference between the nominal angle-of- 
attack profile and the angle-of-attack profile obtained by using the attitude bias from 
figures 4 and 5 is shown in figures 6 and 7. These results were obtained from a detailed 
six-degree-of-freedom computer simulation. The angle-of-attack errors are quite 
small compared to the derived angle-of-attack biases of figures 2 and 3. This indicates 
that the simplified analytic attitude-bias equation is quite accurate. Note that in figures 
6 and 7 the angle-of-attack error using noniterated attitude bias is larger than that for 
the iterated attitude bias. This is caused by the large drift velocities encountered with 
a noniterated bias. The larger the drift velocity, the greater the deviation from the 
nominal trajectory, and the more the accuracy of the equations is reduced. The iterated 
bias anticipates the wind profile. Therefore, the attitude is changed in advance, which, 
in turn, changes the direction of the relative velocity vector. The velocity vector is 
changed in such a way that adding the wind velocity and sideslip effect to it results in 
the nominal relative velocity vector. This minimizes sideslip, which, in turn, improves 
the accuracy of the linearized attitude-bias equation. The noniterated bias does not 
anticipate the wind, and the wind-induced angle of attack must be cancelled by the atti 
tude change. Because the velocity vector is not biased, there is no reduction in side- 
slip, and a larger angle-of-attack error results. 

Figures 8 and 9 give a comparison of altitude dispersions using iterated and non- 
iterated attitude biases. The results were obtained from a six-degree-of-freedom com- 
puter simulation by using the attitude biases derived by the simplified procedure and the 
corresponding wind profile. In all cases, the altitude dispersions obtained by using 
iterated end conditions were much smaller than the altitude dispersions for noniterated 
attitude bias. Figure 8(c) shows a maximum altitude error at 110 seconds of 12 450 
meters for the noniterated bias compared with 150 meters using the iterated bias. This, 
of course, is very significant, since the vehicle is only 36 570 meters in altitude at this 
flight time and a 12 450-meter dispersion can cause severe aerodynamic heating prob- 
lems and payload loss. 

Figures 8 and 9 also give a comparison between a perfect load-relief control sys- 
tem and biased pitch programs. The noniterated attitude bias gives the same flight at- 
titude that a perfect load- relief system would follow. Because of the large dispersions 
introduced, actual load-relief systems are designed to compromise between load reduc- 
tion and flight dispersions. Because biased pitch programs are derived by using the 
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iterated end conditions, the attitude, flight path angle, and attitude dispersions are 
minimized. 

The results presented indicate that the simplified analytic equations give good re- 
sults. In the discussion that follows, the simplified equations are used to compute 
biased pitch and yaw programs. The biased pitch and yaw programs are derived for a 
sample of 100 March wind measurements (ref. 2). The measurements were taken at 
the Eastern Test Range between 1956 and 1959. The pitch and yaw wind components 
were computed based on wind direction and launch azimuth. These wind components 
were used to compute pitch and yaw attitude biases. The launch azimuth for both vehi- 
cle configurations was 105 degrees. 

Figures 10 to 13 give the attitude-bias results for configuration I in the pitch plane. 
The simplified angle- of-attack bias equation was used to obtain angle- of -attack profiles 
for the wind sample. The statistical average of these angle- of-attack profiles is given 
by the dashed curve in figure 10. This angle- of-attack bias was used to derive the at- 
titude bias of figure 11 by applying the simplified attitude-bias equation and the desired 
final condition. To check the accuracy of the bias equation, the attitude bias of figure 11 
was simulated by the six-degree-of-freedom computer program. If a zero wind is used 
in the simulation, the angle- of-attack profile should be equal to the angle- of-attack bias 
of figure 10, except for the inaccuracies of the simplified equations. The resulting 
angle- of-attack bias is given by the solid curve in figure 10. It compares well with the 
statistical angle-of-attack bias. This indicates again that the simplified attitude-bias 
computation gives good accuracy. Since the nominal (unbiased) trajectory is near zero 
angle of attack, the angle-of-attack bias of figure 10 is essentially the same as the 
angle-of-attack history of the biased no- wind trajectory. 

Figure 12 is a histogram of the wind sample analyzed based on maximum angle of 
attack times dynamic pressure. As was expected, the attitude bias drastically reduced 
maximum angle of attack times dynamic pressure, on the averrge. Thus, the probabil- 
ity of having a wind with smaller maximum angle of attack times dynamic pressure is 
increased. This is shown in figure 13. For example, if configuration I had an angle- 
of-attack-times-dynamic-pressure capability of 120 000 degrees times newtons per 
square meter (deg x N/m ) (from fig. 13), 6 percent of the winds considered would not 
exceed this value if the nominal pitch program is used, while approximately 70 percent 
of the winds would not exceed the vehicle capability if the biased pitch program derived 
is used. As another example, if configuration I were designed with 85 percent launch 
availability, its design structural strength would have to provide for a 390 000- 

o 

deg x N/m capability without biased pitch programs. The requirement can be reduced 
c 2 

to 180 000 deg x N/m using the biased pitch program. From the preceding discussion, 
the advantages of using biased pitch and yaw pregrams are obvious. 

Figures 14 to 18, with the exception of figure 16, correspond to figures 10 to 13 but 
for configuration n in the pitch plane. However, this configuration had a biased nominal 
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pitch program. The biased nominal pitch program is given in reference 4, where it is 
designated as the one to be used for the month of March. This pitch program was de- 
rived based on a synthetic wind profile also given in reference 4. Figure 16 shows the 
angle-of-attack profiles for the two biased pitch programs, when the vehicle is flown on 
a no-wind trajectory. Note that the nominal pitch program has a maximum angle of at- 
tack of 3. 4 degrees (fig. 16), while the biased pitch program derived herein has a 
2. 3 degree maximum angle-of-attack bias. Because of this overbias, the nominal pitch 
program gives a reduced launch availability. This is shown in figure 18. For example, 
if configuration II had an angle- of -attack-times-dynamic-pressure capability of 50 000 
deg x N/m (from fig. 18), 18 percent of the winds considered would not exceed this 
value if the nominal pitch program is used, while approximately 52 percent of winds 
would not exceed the vehicle capability if the biased pitch program derived is used. 

To illustrate the applicability of the simplified equations to the yaw plane, a biased 

yaw program was derived for configuration I. The results are given on figures 19 to 22. 

Since the nominal trajectory is zero angle of attack, the angle-of-attack bias is identical 

to the total angle of attack. Therefore, the total-angle-of-attack profile is not plotted. 

o 

From figure 22, if the vehicle capability is 100 000 deg x N/m , the launch availability 
can be improved from 47 percent without bias to 72 percent with the bias. 

To determine the launch availability of a particular vehicle using biased pitch pro- 
grams, the total angle of attack times dynamic pressure must be computed for the wind 
sample. Once the angle- of-attack-times-dynamic-pressure histories are known, the 
same statistical procedure can be used as previously to derive launch availability. The 
total angle of attack can be calculated from the pitch and yaw angles of attack by using 
spherical geometry. 


SUMMARY OF RESULTS 

The results of this analytic calculation of launch vehicle response to winds are two- 
fold: 

1. First, a set of analytic equations was derived to evaluate 

a. The change in angle of attack caused by a wind disturbance for a trimmed 

vehicle (eq. (10)) (A vehicle is trimmed if nominal attitude is maintained 
in the presence of aerodynamic disturbances. ) 

b. The change in attitude caused by a wind disturbance if the nominal angle-of- 

attack profile is maintained (eq. (4)) 

c. The change in attitude necessary to obtain a desired angle-of-attack profile 

(This equation (eq. (5)) was derived by assuming a no-wind nominal trajec- 
tory. ) 
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d. Deflection requirements for a particular angle-of-attack change (eq. (11)) 
The analytic equations can be used in preliminary design studies to determine the 
vehicle's structural strength requirement. This determination can be based on a large 
sample of actual wind soundings instead of a synthetic wind profile. Since the results 
are based on a large number of wind measurements, they will give a more accurate 
representation of load requirements. 

2. Second, the simplified equations were used to derive biased pitch and yaw pro- 
grams. It was shown that substantial improvement in launch availability is obtained by 
using the procedure described herein compared to a biased pitch program derived based 
on an artificial wind profile. The launch availability is computed based on maximum 
angle of attack times dynamic pressure, which assumes the vehicle has a constant angle- 
of-attack-times-dynamic-pressure capability. If the variation of the vehicle capability 
with flight time is known, the launch availability can be computed based on the given 
vehicle capability. Also, the present procedure requires a minimum of computer time 
to derive such a bias. To derive a biased pitch or yaw program using a sample of 
100 winds requires approximately 7 minutes of (IBM 7094) computer time compared with 
approximately 500 minutes using the same wind sample and a six-degree- of-freedom 
computer program. 

The simplified equations may make it feasible to derive a biased pitch and yaw pro- 
gram at the time of launch for the existing launch wind. This has not been evaluated in 
the present report. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 25, 1970, 

731-25. 
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APPENDIX A 


SYMBOLS 

constant defined by eq. (B5b) or (B9b), sec - * 
a 2 constant defined by eq. (B2c) or (B7b), N 

bj constant defined by eq. ,(B5c) or (B9c), sec - * 

normal force coefficient 

C N 0 normal force coefficient for zero angle of attack 
Cj constant defined by eq. (B5d) or (B9d), sec - * 

Ej interval defined by t^ < t < t^ + j 

axial force, N 

Fj interval defined by tj < t < tj + ^ 

F||y force parallel to velocity vector, N 

F^y force perpendicular to velocity vector, N 

o 

g gravitational acceleration, m/sec 

2 

I moments of inertia, N-m-sec 

l „ aerodynamic moment arm, m 

l control moment arm, m 

Mqq moment about center of gravity, N-m 
m gross mass, kg 

N total normal force, N 

Nj normal force per angle of attack, N/rad 

Ng zero-angle- of-attack normal force, N 

O 

Q dynamic pressure, N/m 

Rp Rg defined by eqs. (Bla) and (Bib), N 

2 

S re f aerodynamic reference area, m 
T total thrust, N 

t time, sec 

V velocity, m/sec 

W g effective gravity force, N 



x time derivative of x, sec" 

a angle of attack, rad 

y flight path angle, rad 

A linearized variable (i. e. , AX = X - X n ) 
6 thrust deflection angle, rad 

9 flight attitude, rad 

9 reference pitch attitude, rad 

Jr 

£ dummy variable 

r dummy variable 

Subscripts: 
e effective 

f final 

i, j time intervals 

n nominal 

o starting value 

r relative 

w wind 
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APPENDIX B 


DERIVATION OF BASIC EQUATIONS 

The basic equations used in the analysis will be derived for a rigid-body vehicle 
configuration. The pitch, yaw, and roll planes are assumed to be uncoupled; that is, 
disturbances in any of these planes will not be felt in the other two planes. This is a 
good approximation for symmetric or nearly symmetric vehicles. The autopilot effects 
(the effects of a real control system) are not included in the equations. Equations are 
derived separately for the pitch and yaw planes. However, the final results are very 
similar, and the same equations can be used with minor modifications to analyze both 
planes. The variables used are defined in appendix A. 


Equations in Pitch Plane 

Since the vehicle is roll- stabilized, the pitch plane is defined by the radius vector 
and the local horizontal in the launch azimuth direction. The basic vehicle configuration 
in the pitch plane is given in figure 23. The equations of motion are 

^F ||V = Rj = mV = T cos(0 - y + 6) - cos(0 - y) - N sin(0 - y) - mg Sin y 

(Bla) 

^F 1V = R 2 = mVy = T sin(0 - y + 6) - F^ sin(0 - y) + N cos(0 - y) - mg cos y 

(Bib) 

£m cg = 10 = NJ a - T£ c sin 6 (Blc) 

a = 0 - y - oiyj (Bid) 

V 

sin cl.. = — — sin y (Ble) 

V 

r 

where V w is the wind velocity component in the pitch plane, and 

N = Njo; + Ng (Blf) 

where the normal force per angle of attack = QS re j(dCj^/daf), and the normal force 
due to the asymmetry of the vehicle Ng = q. 
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Linearizing equation (Bib) and evaluating the coefficients along the nominal trajec- 
tory (which is assumed to be a no-wind trajectory) gives 

my n AV + niV n Ay = ARg (B2a) 


where 


ARg = ^0 ~ (R]) A y + cos « n A O' + T cos(o! + 6) n A6 (B2b) 

a .2 = T cos (a + 6) n - cos <* n - N sin a n (B2c) 

The nominal case is assumed to be a no- wind trajectory. Since the change in relative 
velocity from the nominal caused by the wind disturbance is small, the term my n AV in 
equation (B2a) can be neglected compared to mV n Ay, and the equation reduces to 

AR 9 

Ay = “ (B3) 

mV n 


Linearizing equations (Blc) and (Bid) gives 

A« = ! (N.l Aa - I A0) (B4a) 

K c c 0S6 n 

Ay = A0 - A a - (B4b) 

Eliminating A5 and Ay in equations (B2b), (B3), (B4a), and (B4b) gives the change in 
attitude in terms of the change in angle of attack and wind- induced angle of attack; that 
is, 


I cos( oi + 6) 


m Vc cos 6 n 


n ** 

— A0 + A0 + a^ A0 =Aa + c 1 Aa + + b l 


% 


This equation can be reduced further by assuming the inertial effects of the vehicle on 
the change in attitude to be small (i. e. , 1 = 0). Then, 


A £ = a^ A 9 = A a + Cj Aa+ + b^a w 


(B5a) 
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where 


W. 


a l = - 


mV 


n 


(B5b) 


(Ri) 


v 1 ; 

b i =— B = c °s(« + «) n " F A cos - N sin « n - wj 


mV mV 
n n 


(B5c) 


c l = b l + 


Ni 


mV 


n 


cos + 


Z a cos(a + 6) n 

l n cos 5_ 
c n _ 


(B5d) 


W e = mg sin y fi 


(B5e) 


Equations (B5) are the basic equations used in the analysis. 


Equations in Yaw Plane 


The yaw plane is defined by the vector perpendicular to the launch azimuth plane at 
the vehicle’s center of gravity and the projection of the longitudinal axis in the azimuth 
plane. The basic vehicle configuration in the yaw plane is given on figure 24. The 
equations of motion are 


F,| V = Rj = mV = T cos(0 - y + 5) - F A cos(0 - y) - N sin(0 - y) - mg sin 6 p cos y 

(B6a) 

F iV = r 2 = mVy = T sin(0 - y + 5) - F^ sin(0 - y) + N cos(0 - y) + mg sin 0 p sin y 

(B6b) 


M CG = B=N 'a- T 'c sin5 


oi — 9 - y - of, 


w 


sin a_ = — — cos y 

W y 

r 


(B6c) 

(B6d) 

(B6e) 
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N = N 1 q ! + N 2 


(B6f) 


Linearizing equation (B6b) and assuming the change in velocity caused by wind disturb- 
ance to be small gives 

Ay = -1— AR 9 (B7a) 

mV n 


where 


ARg = a 2 A0 - (Rj) Ay + cos A or + T cos(oi + S) n A6 (B7b) 


a 2 = T cos( a + 6) n - cos cn n - N sin a n (B7c) 


Linearizing equations (B6c) and (B6d) gives 


A6 = - (N,Z Aa - I A0) (B8a) 

t; c cos6„ 


Ay = A0 - A a - ot^ since (a^) = 0 (B8b) 

n 

Combining equations (B7) and (B8) and assuming the inertial effect to be negligible gives 
the attitude change in terms of the change in angle of attack and wind-induced angle of 
attack 

• • • 

A 0 + a^ A0 = Aof + Cj Aa; + + bjQ! w (B9a) 


where 


W e 

a, = 5- (B9b) 

mV n 


(Rj) 

b-i = = — — [t cos(a + 6)„ - F A cos a - N sin ol. - W_1 (B9c) 

1 mV mV. L n A n n eJ 

n n 
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W g = mg sin 0 p cos y n (B9e) 

Note that these equations are identical with equations (B5) derived in the pitch plane, 
with the exception of the equivalent weights given by equations (B5e) and (B9e). There- 
fore, the same equations may be used in the analysis in both pitch and yaw planes by 
using the appropriate equivalent weights given by equations (B5e) and (B9e), respec- 
tively. 
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Figure 3. -Change in angle of attack from nominal, for configuration II. 
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(a) Using wind 1. 
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Figure 4. - Simplified analytic attitude bias from nominal, for configuration I. 
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(b) Using wind 2. 

*e 5. - Simplified analytic attitude bias from nominal, for configuration II. 
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Figure 10. - Pitch plane angle-of-attack bias obtained by using the simplified analytic equations and 
the sample of 100 March winds, for configuration I. 
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Figure 11. - Simplified analytic attitude bias using angle-of-attack bias of figure 10, for configuration I. 
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Launch availability, percent 



Maximum angle of attack times dynamic pressure, deg x N/m' 


Figure 12. - Histogram for distribution of 100 March winds obtained by using simplified ana- 
lytic equations, for configuration I. 



Maximum angle of attack times dynamic pressure, deg x N/nr 


Figure 13. - Launch availability derived by using simplified analytic equations with and without attitude bias of fig- 
ure 11, for configuration I. 
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Flight time, sec 

Figure 14. - Pitch-plane angle -of -attack bias obtained by using the simplified analytic bias and the 
sample of 100 March winds, for configuration II. 
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Figure 15. - Simplified analytic attitude bias using angle-of-attack bias of figure 14, for configuration I 





Figure 19. - Yaw plane angle-of-attack bias obtained by using the simplified analytic equations and the 
sample of 100 March winds, for configuration I. 
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Figure 20. - Simplified analytic attitude bias in yaw plane using angle-of-attack bias of figure 19, for 
configuration I. 



Figure 21. - Histogram for distribution of 100 March winds obtained by using the simplified analytic equa- 
tions, for configuration I. 
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Figure 22. - Launch availability derived by using simplified analytic equations with and without attitude 
bias of figure 20, for configuration I. 
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